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ABSTRACT 

Recent experiments suggest that the heat-transfer 
characteristics of stator blades are influenced by the 
frequency of passing of upstream rotor blades. The 
calculation of these effects requires that the movement 
of the stagnation point with variations in freestream 
velocity is properly represented together with the 
possible effects of turbulence characteristics on the 
thin leading-edge boundary layer. A procedure to 
permit the achievement of these purposes is described 
for laminar flows in this paper together with results 
of its application to two model problems which demon- 
strate its abilities and quantify the influence of wake 
characteristics on fluid-dynamic and heat-transfer prop- 
erties of the flow and their effects on surface heat 
transfer. 


w 

x,y 

x,y 

z,z Q 

a 

Y 

A* 

n 

v 

l 

P 

T 

4> 

* 


= Ug/Ua, 

= boundary- layer coordinates 
= Cartesian coordinates 
= 4>/y f a/y 
= angle of attack 
= thickness ratio = b/a 

- dimensionless displacement thickness, (17b) 
= (u 0 /vL) 1/2 y 

* kinematic viscosity 
= X/L 
= density 

= nondimens ional time, tu 0 /L 
= polar angle related to x,y coordinates 

- stream function 
= frequency 


Subscripts and Symbols 


NOMENCLATURE 

a,b = ellipse axes 

D = cylinder diameter 

F = blade passing frequency, l/(t g + t w ) 
f = stream function, i|r/(u 0 Lv) 

G = nondimens ional temperature, (T w - T)/(T W - T e ) 
k = thermal conductivity 

L = reference length 

m = (T w - T e ) -1 [d(T w - T e )/d£] 

Pr = Prandtl number 

q^ = wall heat flux 

R™ * Reynolds number, u c L/v 

St = Strouhal number, FD/uo, 

T = static temperature 

t = time 

u,v = boundary- layer velocity components 

u e = u e /[u«»(l + y) ] 

u Q = reference velocity, equal to u«,(l + y) and U® 
for first and second model problems, 
respectively 

u® = freestream velocity 


e = boundary layer edge 

s = stagnation point 

w = wall or wake 

primes denote differentiation with respect to n 
INTRODUCTION 

Periodic reviews of gas turbine heat transfer 
( 1-3 ) suggest that progress is being made in measuring 
and analyzing the complex flows associated with steady- 
state heat transfer in gas turbines. However, research 
on unsteady heat transfer in turbine passages is just 
beginning to appear. It is premature to evaluate the 
effects of wake-generated unsteady heat transfer on 
turbine durability but the research data to provide the 
necessary thermomechanical loads ( 4-13 ) are emerging 
rapidly. Dring, et al. (4) explored the nature of the 
boundary-layer response to turbine wakes and the un- 
steady pressure loading using thin-film surface sensors. 
Hodson (6) and Binder, et al. (7) mapped the convection 
of wakes through the rotor and flow properties in a 
rotor. Most of the heat-transfer data have been 



obtained in short-duration facilities ( 8-11 ) or by using 
short duration techniques (12,13). Thin-film resistance 
thermometers provide information of surface temperature 
and are combined with a one-dimensional semi-infinite 
heat-conduction analysis to yield time-resolved heat 
flux. Dunn and coworkers have pioneered this technique 
and, in their most recent work (8,9), have begun to 
obtain time-resolved heat flux data in a real turbine 
stage at fully-scaled conditions. At the University of 
Oxford (10*11) rotor wakes are simulated with a spoked 
wheel rotating in front of a linear cascade. The re- 
sults show that the wakes impinging on the airfoil can 
cause the local boundary layer to undergo transition and 
reverse transition at high frequency of the passing 
wakes. O'Brien et al. (12) used the transient technique 
in a steady-running rig, also with a spoked wheel wake 
generator, and showed excellent agreement between the 
time-average of the time-resolved heat transfer and more 
conventional steady-state data. Morehouse and Simoneau 
( 13) in the same rig showed that it was necessary to 
separate the heat transfer effects due to the unsteady 
disturbances of periodic wakes from those generated by 
conventional turbulence. 

In order to characterize the complex time- 
dependent flowfields, a number of promising approaches 
are being developed. Hodson (14) employed the two- 
dimensional inviscid code of Denton (15) and imposed a 
quasi-steady inlet boundary condition of a wake shape 
which moved across the inlet of the cascade at blade 
speed and calculated the movement of the wake as it 
passed through the cascade. Rai (3JO has constructed a 
Navier-Stokes computer code with a complex arrangement 
of patched and overlaid grids and has attempted to re- 
produce the measurements of Dring (4) with some success. 

An alternative approach which has proven reliable 
and very powerful in many external flows is to use 
interactive boundary- layer theory which requires the 
development of appropriate viscous and inviscid flow 
methods and their coupling by special techniques such 
as those described in ( 17-21) . 

As an essential preliminary to more extensive 
calculations involving complete blades, here we consider 
the development of an unsteady boundary-layer method for 
calculating the flow properties near the stagnation 
region of a blade where the movement of the stagnation 
point with space and time poses problems. We assume 
that the external velocity distribution is represented 
in dimensionless form by a function 


BASIC EQUATIONS 


For a two-dimensional, incompressible time- 


dependent laminar flow, the boundary-layer equations 
and their boundary conditions are well known, and for 
conditions with no mass transfer and specified wall 
temperature they can be written in the form: 


1“ + f* = 0 

3x 3y 


(2) 
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ST ^ 3T 3T V 3 2 T 
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(4) 


y = 0, u = v = 0, T = T (x) ; 
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(5a) 


y = 6, u=u e (x,t), T=T^ 


(5b) 


The determination of the initial conditions 
required for the above system is important and sometimes 
it can be arbitrary but in that event , the values of 
3u/3t at t = 0 are nonzero; this implies an invis- 
cid acceleration and, as a consequence, a slip velocity 
develops at the wall and is smoothed by an inner bound- 
ary layer initially of thickness (vt) 1/2 in which 
viscous forces are important. Thus a double structure 
develops in the boundary layer and may be treated by the 
numerical method described in (22). However, if inter- 
est is centered on the solution at large times, this 
feature may be reduced in importance by requiring that 
the initial velocity distribution satisfies the steady- 
state equation with the instantaneous external velocity. 
In addition, it is necessary to smooth out the external 
velocity u e (x,t) so that 3u e /3t = 0 at t = 0 and then 
standard numerical methods may be used and are stable. 
The use of a smoothing function makes for some loss of 
accuracy at small values of t but the error soon decays 
to zero once the required value of u e is specified. 


The calculation of upstream boundary conditions in 
the (t,y) plane at some x = x Q when the conditions at 
a previous time line are known, can introduce different 
problems. To illustrate these difficulties for the case 
of a moving stagnation point, let us consider Eq. (1). 
Since u e = 0 at the stagnation point by definition, 
its location based on the external streamlines is given 
by 

l s = B(t) (6) 


U e /Ua, (X,t) = A(£,l)[£ - B(t) ] (1) 

which allows the variation of the stagnation point and 
the freestream velocity. To evaluate the approach and 
the numerical procedure, we have chosen two model prob- 
lems. The first one corresponds to flow over a thin 
ellipse in a constant freestream velocity field with a 
uniform wall temperature and with a local external vel- 
ocity which changes with angle of attack. The second 
model problem considers flow over a circular cylinder 
again with a uniform wall temperature but with an exter- 
nal velocity distribution which allows changes in free- 
stream velocity with time and simultaneous movement of 
the stagnation point caused by the change in position 
of a wake from a rotor blade. 

The basic equations, initial and boundary condi- 
tions are considered in the following section and are 
followed by the results section. The paper ends with a 
summary of the more important conclusions. 


Figure 1 shows the variation of the stagnation 
point with time according to the above equation with 
B(t ) = 1 + c sinwT, c = 1, G) = ir/4. We see that the 
stagnation point l s is at 2 when t = -2 and at 0 when 
t = 6, etc. If l s were fixed, we could assume that u = 0 
at l s = -1 for all time and for all y, but this is 
not the case. It is also possible to assume that the 
stagnation point is coincident with zero u-velocity for 
a prescribed time. However, we should note that the 
stagnation point given by Eq. (6) is based on vanishing 
external velocity. For a time-dependent flow, as we 
shall discuss later, this does not necessarily imply 
that the u-velocity is zero across the layer for a given 
^-location and specified time; flow reversals do occur 
due to the movement of the stagnation point and cause 
the locus of zero u-velocity to vary with l requiring 
the use of a special numerical method. 

It is more convenient and useful to express Eqs. 
(2)-(5) in a form more suitable for computation. To 
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Figure 1. Variation of the stagnation point with time 
according to Eq. (1). 


achieve this, we introduce the dimensionless variables 
t, l, Hr w, ra, G together with a dimensionless stream 
function f(£,n»T) and, with 9 = 3f/3£, write Eqs. (2) to 
(5) as 


+ f 


■ n 3w , dw _ 

6 + 3t + w 35 " 3 t 


+ f' 


3f ' 

35 


(7) 


and Pr = 0.72. In regions of no flow reversal, the 
Standard Box scheme was used. Where a calculation with 
this scheme revealed flow reversal, further iterations 
at that location made use of the Characteristic Box 
scheme described in (24-26 ) for two-dimensional unsteady 
flows and in (23,27) for three-dimensional steady flows. 
The application of the Characteristic method to the 
present problem, together with the procedure used to 
generate the upstream conditions is described in (28) 
and for this reason only the results are presented and 
discussed here for the two model problems. In all the 
calculations done here, the initial conditions were 
obtained by solving the steady-state form of the equa- 
tions with the Standard Box scheme for uniform step- 
lengths in l and in n. The stagnation point solu- 
tions were obtained from the similarity equations. A 
novel procedure based on the Characteristic Box scheme 
was used to compute the first solution on each time line 
with the initial calculations starting near the stagna- 
tion point based on the external velocity given by Eq. 
( 1 ). 


— G" + G * 0 + m( 1 - G)f * - + f* If (8) 

Pr ox 

n - 0, f = f* = G = 0; n = n e , f* = w, G = 1 (g) 

We use the Characteristic Box method to avoid the 

numerical problems associated with flow reversals in the 
stagnation region and to generate the initial conditions 
on the next time "line." This scheme requires that Eqs. 
(7) and (8) are expressed in terms of new coordinates. 
For this purpose, we note the definition of local 

streamlines and let di = d£/f'. If the distance in 
this direction is designated by s and the angle that it 
makes with the t-axis by P, (see Fig. 2) then the 

transformed momentum and energy equations (7) and (8) 
and their boundary conditions can be written as 


Oscillating Airfoil Model 

The first model problem corresponds to flow over 
an ellipse with a thickness ratio y(Eb/a) and with 
Y « 1 at an angle of a. The surface of the body is 
defined by 

x = -a cos4>, y = ay sin<|>, -if < 4> < ir (14) 

With these definitions and to a first-order approxima- 
tion, the external velocity for the steady flow in the 
leading edge region of a thin ellipse can be deduced 
from inviscid flow theory to be 

z + z 

U (z) = — — — (15) 

/ 1 + z 2 


1_ 

Pr 


where 


* '■* * If * » if ■ *- !r 

(10) 

+ G’9 + m(l - G) f * = X 

(11) 

= f' = G = 0, n = n , f = w, 
e 

G = 1 (12) 

+ (f') 2 , P = tan _1 f ' 

(13) 


RESULTS AND DISCUSSION 


Here u G (z) denotes a dimensionless velocity, the 
parameter z _denotes a dimensionless distance related to 
the x- and y-coordinates of the ellipse by x + a = 
(l/2)ay 2 z 2 , y - ay 2 z measured from the nose, and z Q rep- 
resents a reduced angle of attack. The parameter z is 
also related to the dimcris ion less surface uioLan^c £ 
with L = y 2 a by l = oJ 2 (l + z 2 )^^ 2 dz. in the present 
study we extend Eq. (15) to unsteady flows by introduc- 
ing time dependency as 


u (z,t) 
e 


z + z (1 + c sinurr) 

o 


y i + z 2 


(16) 


The numerical solution of Eqs. (10) and (12) is 
obtained with Keller's Box method discussed, for 
example, in (23) for m = 0 (uniform wall temperature) 



Figure 2. Notation and grid for the Characteristic Box 
scheme. For details, see (28). 


Figures 3 to 5 show the calculated results for two 
circular frequencies, w = ir/30, ir/3 with z Q = 1 
and c = -1/2. Figures 3 through 5 show the effects of 
the frequency on the wall shear f"(0), wall heat flux 
G ' ( 0 ) and dimensionless displacement thickness, A*, 
which are defined by 

t qL 

f"( 0) = -^ VR , G ' ( 0 ) - — (17a) 

pu Q k^(T w - T e ) 

- l ir s (1 - r )dy (17b) 

o 0 e 

The results of Figures 3a and 3b indicate that the 
time dependent effects on the solutions increase signif- 
icantly as the circular frequency cj changes from u =ir/30 
to ir/3. While the wall shear and displacement thickness 
values computed at un = t and 2tr with w = ir/30 are 
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nearly the same as those at x = 0, they differ con- 
siderably from each other when u) = ir/3. The results in 
Figure 4 show, however, that increasing the circular 
frequency by a factor of ten has practically no effect 
on the wall heat flux. 

Figures 5 and 6 allow the examination of the effect 
the frequency has on the calculated velocity profiles 
in the vicinity of the stagnation point. Figure 5 shows 



Figure 3. The effect of frequency on the (a) wall shear 
parameter, f"(0), (b) displacement thickness, 
A*, of the first model problem. 
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Figure 5. variation of the velocity profiles for w = 

ir/30 near the stagnation region of the first 
model problem for different values of wx, (a) 
90, (b) 180, (c) 270, (d) 360. 
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Figure 4. Variation of the wall heat flux parameter, 

G * ( 0 ) , with z for the first model problem 
with a) - tt/ 30. The results with u> = ir/3 
are virtually identical to those obtained 
with the smaller frequency. 


Figure 6. Variation of the velocity profiles for w = 
ir/3 near the stagnation region of the first 
model problem for different values of cjt , (a) 
90, (b) 180, (c) 270, (d) 360. The dashed 

line indicates the locus of zero u-velocity. 

that the locus of the u-velocity on time lines x = ir/2 
and ir is essentially the same as in the steady case and 
as a result there are no flow reversals in the velocity 
profiles. However, as can be seen from Figure 6, in- 
creasing the frequency to ir/3, flow reversals begin 
to occur around the stagnation point and become rather 
prolonged as time increases to wx = 3ir/2. At cox = 2ir, 
the region of flow reversal is reduced but is not zero, 
as it was at ux = 0. 
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Rotor Wake Model 


The second model problem corresponds to a flow on a 
circular cylinder of diameter D, which experiences the 
periodic passing of the wakes from the turbine blades, 
see Fig. 7. For a time period t g , the cylinder is 
subjected to a freestream velocity and for t w it is 

immersed in a superimposed moving wake which has a 
motion component u^. The cycle repeats itself with a 
so-called blade passing frequency F[i l/(t g +t w )] 
and is related to the strouhal number St by FD/u«,. 
We assume that the stagnation region of the cylinder is 
subjected to a velocity which varies in space and time 
according to Eq. (1), see (28), where 

R(T) = [h 2 + 10(St) 2 (l - h) 2 ] 1/2 

B(t ) = tan -1 [3.3St (i - X)] 

1 + u 1 - U 

h = - 2 m + 2 m cos[tr(l - 100 s t t ) } 

Here, u m denotes the dimensionless minimum 
velocity. 


(18) 

(19) 

(20) 
wake 


Calculations are made for two values of strouhal 
number taken equal to 0.1 and 0.2 with 1 % * 1/3 in 
both cases. The computed values of wall heat flux, 
G’(0) show that they are not influenced by the changes 
in the freestream velocity and are virtually constant 
for the range of £ and t values considered with 
G'(0) - 0.50 for St = 0.1 and G*(0) ~ 0.51 for St = 0.2. 



Figure 7. Notation and flow configuration for the 
second model problem. 


On the other hand, as shown in Fig. 8, the computed 

values of wall shear, f*(0), for St - 0.1 are signifi- 
cantly influenced by the changes in the freestream vel- 
ocity which causes flow reversals in the velocity pro- 
files around the stagnation point based on the vanishing 
of the external velocity. The movement of the stagna- 
tion point and the resulting flow reversals increase 
with time and with space. For example, the calculations 
for steady state have the stagnation point at £ * 0, 
and, as expected, there is no flow reversal in either 
side of the stagnation point. At t = 0.05, the stag- 
nation point moves to £ = 0.15 but the flow reversals 
in the velocity profiles continue up to and including 

f - 0.85 as can be seen from the results shown in 

Fig. 9. At t = 0.10, Fig. 9, the stagnation point 

has moved to £ = 0.55 but the flow reversals persist 

for a longer distance and continue until \ = 1.20. 

As can be seen from the velocity profiles in Fig. 9, 

the region of flow reversal across the layer has now 
increased and is substantially more pronounced than 

those at t = 0.05. 



Figure 8. Variation of the wall shear parameter, f"(0), 
with \ for the second model problem. St = 
0 . 10 . 




Figure 9. Variation of the velocity profiles near the stagnation region of the second model problem for two values 
of time, (a) t = 0.05, (b) t = 0.10. 
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The wall shear and displacement thickness results 
for St = 0.2 show similar trends. As expected, the 
flow reversals in the velocity profiles, for example at 
t = 0.10, are bigger than those at St - 0.1 but cover 
the same range in £. 

CONCLUDING REMARKS 


A method for computing the heat transfer in the 
stagnation region of turbine blades is described for 
incompressible flows and evaluated for two model prob- 
lems involving laminar flows. This method unlike the 
similarity methods of (29-31), however, is general and, 
with minor changes, can also be used for turbulent flows 
as well as for compressible flows. 

The computed results show that the numerical pro- 
cedure is able to obtain solutions for a range of blade- 
passing frequencies of practical relevance. The move- 
ment of the stagnation point with space and time and the 
resulting flow reversals around the stagnation point 
cause no computational difficulties and the numerical 
tests show that the accuracy is better than required for 
practical problems. 

Acknowledgment: The research reported in this paper was 
conducted under NASA Lewis Research Center Grant, MIPR 
C-80017-F. 
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